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Abstract 

A detailed description of the CUORICINO 130 Te neutrinoless double-beta (Oiy/3/3) decay experiment 
is given and recent results are reported. CUORICINO is an array of 62 tellurium oxide (TeC>2) bolometers 
with an active mass of 40.7 kg. It is cooled to ~ 8 — 10 mK by a dilution refrigerator shielded from 
environmental radioactivity and energetic neutrons. It is running in the Laboratori Nazionali del Gran 
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Sasso (LNGS) in Assergi, Italy. These data represent an exposure of 11.83 kg • y or 91 mole-years of 
130 Te. No evidence for O/3/3-decay was observed and a limit of T°J 2 ( 130 Te) > 3.0 x 10 24 y (90% 
C.L.) is set. This corresponds to an upper limit on the effective mass, (m„), between 0.19 and 0.68 eV 
when analyzed with the many published nuclear structure calculations. In the context of these nuclear 
models, the values fall within the range corresponding to the claim of evidence of O^/3/3-decay by H.V. 
Klapdor-Kleingrothaus, et al. The experiment continues to acquire data. 



I INTRODUCTION 

There are three very important open questions in neutrino physics that can best be addressed by next gener- 
ation neutrinoless double-beta (0v/3(3) decay experiments. First, are neutrinos Majorana particles that differ 
from antineutrinos only by helicity? Second, what is their mass-scale? Third, is lepton number conservation 
violated? While searches for /3/3-decay have been carried out steadily throughout many decades [1-3], it is 
now a far more interesting time for the field. Atmospheric neutrino-oscillation data imply that there exist 
scenarios in which the effective Majorana mass of the electron neutrino could be larger than 0.05 eV. Recent 
developments in detector technology make the observation of O^/3/3-decay at this scale now feasible. For 
recent comprehensive experimental and theoretical reviews see [4-6]. Optimism that a direct observation of 
Oz^/3/3-decay is possible was greatly enhanced by the observation and measurement of the oscillations of at- 
mospheric neutrinos [7], the confirmation by SuperKamiokande [8] of the deficit of 8 B neutrinos observed 
by the chlorine experiment [9], the observed deficit of p — p neutrinos by SAGE [10] and GALEX [11], 
and the results of the SNO experiment [12] that clearly showed that the total flux of 8 B neutrinos from the 
sun predicted by Bahcall and his co-workers [13] is correct. Finally, the data from the KamLAND reactor- 
neutrino experiment strongly favor the MSW large mixing-angle solution of solar neutrino oscillations [14]. 
This important list of results published since 1998 weighs very heavily in favor of supporting two or more 
next generation O^/3/3-decay experiments (see the reports in references [15, 16]). 

The most sensitive limits have come from germanium detectors enriched in 76 Ge. They were 
the Heidelberg-Moscow experiment (t"J 2 ( 76 Ge) > 1.9 x 10 25 y) [17] and the IGEX experiment 

(^1/2 ( 76( ^ e ) ^ 1-6 x 10 25 y^ [18]. These imply that the upper bound on the effective Majorana mass 
of the electron neutrino, (m„), defined below, ranges from ~ 0.3 to ~ 1.0 eV, depending on the choice of 
nuclear matrix elements used in the analysis. However, a subset of the Heidelberg-Moscow Collaboration 
has reanalyzed the data and claimed evidence of a peak at the total decay energy, 2039 keV, implying Quf}f3- 
decay [19,20]. While there have been opposing views [21-23], there is no clear proof that the observed peak 
is not an indication of O^/3/3-decay. The GERDA experiment, also using 76 Ge, is under construction in the 
Laboratori Nazionali del Gran Sasso (LNGS), and will test this claim [24]. The CUORICINO experiment, 
also located at LNGS, is the most sensitive O^/3/3-decay experiment with good energy resolution currently 
operating [25,26]. It is searching for the O^/3/3-decay of 130 Te and has the capability of confirming the 
claim; however, a null result cannot be used to refute the claim because of the uncertainty in the nuclear ma- 
trix element calculations. The proposed Majorana 76 Ge experiment [27], CUORE 130 Te experiment [28], 
and EXO 136 Xe experiment [29] are all designed to reach the (m u ) 0.05 eV mass sensitivity and below. 
Descriptions of other proposed experiments with similar goals are given in the recent reviews [4-6]. 

There are other constraints on the neutrino-mass scale, irrespective of their Majorana or Dirac character. 
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The Troitsk [30] and Mainz [31] 3 H single /3-decay experiments have placed an upper limit of 2.2 eV on 
the mass of the electron neutrino. The KATRIN experiment, a greatly enlarged 3 H /3-decay experiment in 
preparation, is projected to have a sensitivity of 0.2 eV [32]. 

Astrophysical data are also very relevant in a discussion of neutrino mass. In a recent paper by Barger et 
al, [33] an upper limit on the sum of neutrino mass eigenvalues, S = mi + + < 0.75 eV (95% 
C.L.), was derived. The data used were from the Sloan Digital Sky Survey (SDSS) [34], the two degree Field 
Galaxy Red Shift Survey (2dFGRS) [35], and the Wilkinson Microwave Anisotropy Probe (WMAP) [36], 
as well as other CMB experiments and data from the Hubble Space Telescope. Hannestad [37] used the 
WMAP and 2dFGRS data to derive the bound E < 1.0 eV (95% C.L.) and concluded that these data alone 
could not rule out the evidence claimed in [19,20]. On the other hand, Allen, Schmidt and Briddle [38] 
found a preference for a non-zero neutrino mass, i.e., S = 0.561q25 This is interestingly close to the 
favored range of values given in [19, 20]. For recent papers on the subject see [39] and references therein. 
The constraint E < 0.75 eV would imply that the lightest neutrino eigenstate mass mi < 0.25 eV. On the 
other hand, if the claim of the positive value of E would be correct, (m„) w 0.17 eV, and next generation 
O^/3/3-decay experiments would constitute a stringent test of lepton-number conservation, irrespective of the 
neutrino mass hierarchy (see the discussion of hierarchy below). 

In this paper we present a detailed description and present the results from the CUORICINO O^/3/3-decay 
experiment derived from data taken between April 2003 and May 2006. Finally, we note that 130 Te has 
a series of calculated matrix elements implying values of (m„) derived from the CUORICINO half-life 
limit between ~ 0.20 keV, and ~ 0.68 keV. A detailed discussion of the implications from the recent 
developments in the theoretical nuclear structure calculations is given later. 

II NEUTRINO PHYSICS AND NEUTRINOLESS 
DOUBLE-BETA DECAY 

Neutrino-oscillation data very strongly imply that there are three neutrino flavor eigenstates, \v e ^_ T ), that 
are super positions of three mass eigenstates, \v\ 2 3), of the weak Hamiltonian as expressed in equation 
Q: 

3 

l^> = El u «l e ^N' (1) 

3=1 

where I = e, fx, r, and the factor e %5 i is a CP phase, ±1 for CP conservation. 

The decay rate for the 0^/3 /3-decay mode driven by the exchange of a massive Majorana neutrino is expressed 
in the following approximation: 

(r^y 1 = G 0v (E ,Z) 



Mf-{g A /gvyM\ 
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where G 0v is a phase space factor including the couplings, |(m„)| is the effective Majorana mass of the 
electron neutrino discussed below, M9 V and M^ T are the Fermi and Gamow-Teller nuclear matrix elements 
respectively, and 17,4 and gv are the relative axial-vector and vector weak coupling constants respectively. 
After multiplication by a diagonal matrix of Majorana phases, (m„) is expressed in terms of the first row of 
the 3x3 matrix of equation (Q~|) as follows: 



l(m,> I 



mi 



m,2e 



i<t>2 



l e3) 



»(<fe+<5) 



(3) 



where e 1 ^ 2 - 3 are the Majorana CP phases (±1 for CP conservation in the lepton sector). Only the phase 
angle 6 appears in oscillation expressions. The two Majorana phases, e 1 * 2 - 3 , do not, and hence do not 
affect neutrino oscillation measurements. The oscillation experiments have, however, constrained the mixing 
angles and thereby the coefficients ufj in equation ([3J. Using the best-fit values from the SNO and Super 
Kamiokande solar neutrino experiments and the CHOOZ [40], Palo Verde [41] and KamLAND [14] reactor 
neutrino experiments, we arrive at the following expression in the case of the normal hierarchy: 



|(m„)| = |(0.70+°;°2) mi + (0.30+°;^) m 2 e* + (< 0.05) m 3 e^+ 5 ) 



(4) 



where the errors are approximated from the published confidence levels (C.L.). The bound on |xt e 3 1 is at 
the 2<T C.L. and the errors on the first two coefficients are la. In the convention used here, the expression 
for the inverted hierarchy, discussed below, is obtained by exchanging mi 7713 in equation (0|. 



The results of the solar neutrino and atmospheric neutrino experiments yield the mass square differences 
7if — m|| but cannot distinguish between two mass patterns (hierarchies): the "normal" hierarchy, 



8?j = \mj-m 2 j 
in which Sm^ 



mf and mi = 7712 -C ma, and the "inverted" hierarchy where Sm 2 oW = m 



and ma = to 2 ^> mi. In both cases we can approximate 8m 2 AT = m 2 — m\ . Considering the values in 
equation (0|, we make the simplifying approximation (u e s) 2 « 0. Using the central values of equation (@), 
we can write the following approximate expressions: 



Km,,)! = mi 



0.7 + 0.3e 




(5) 



for the case of "normal" hierarchy, and, 



|(m„)| - + 5m 



2 

AT 



0.7 + 0.3e 



102 



(6) 



in the "inverted" hierarchy case. At this time there is no experimental evidence favoring either hierarchy. In 
Table U we use Eqs. © and © to show the predicted central values of (m„) as a function of the lightest 
neutrino mass eigenvalue, mi. These values roughly define the desired target sensitivities of next generation 
0z//3/3-decay experiments. 
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Table I: Central values of the numerical predictions of | (m„) | (meV) for both hierarchies and CP phase 
relations. ( mi is also given in meV.) 



Normal Hierarchy Inverted Hierarchy 
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24.0 
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32.0 
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80.0 


75.0 


36.0 


75.0 


90.1 


100.0 


40.0 


100.0 


100.0 


100.0 


44.7 


100.0 


111.0 


200.0 


80.0 


200.0 


200.0 


200.0 


82.5 


200.0 


206.0 


400.0 


160.0 


400.0 


400.0 


400.0 


161.1 


400.0 


403.0 



It is clear that a next generation experiment should have at least the sensitivity for discovery in the case of 
an inverted hierarchy when e 1 ^ 2 = e 1 ^ 3 and for mi = 0. In this case, (m^) i=s J <^ T w 0.050 eV. It should 
also be capable of being expanded in case this level is reached and no effect is found [15, 16]. 

It is convenient to define the nuclear structure factor, Fn, (sometimes denoted as C mm in the literature) as 
follows: 



Fjy = G" 



Mf~{g A /g v YM, 



2 „ /r 0i/ 
GT 



(7) 



Accordingly, the effective Majorana mass of the electron neutrino is connected to the half-life as shown in 
equation (|8j: 



(m„) = ; e ■ (8) 

To extract values of F V from theoretical papers, we recommend using their calculated values of half lives 
for a given value of (m v ), thereby avoiding difficulties associated with conventions used in calculating 
phase-space factors. 

Possible interpretations of the null result of CUORICINO, in terms of the effective Majorana neutrino mass, 
may be understood with detailed analyses of the nuclear matrix elements discussed in a Sees. IVIIII andllXI 
In Sec. |X] this null result will be compared with the positive claim report reported in [19, 20]. 
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Ill THE EXPERIMENT 



The CUORICINO experiment is an array of cryogenic bolometers containing Te, the parent O^/3/3-decay 
isotope. This technique was suggested for /3/3-decay searches by Fiorini and Niinikoski [42] and applied 
earlier by the Milano group in the MIBETA experiment [43]. The bolometers are sensitive calorimeters 
that measure the energy deposited by particle or photon interactions by measuring the corresponding rise in 
temperature. The CUORICINO bolometers are single crystals of Te0 2 ; they are dielectric and diamagnetic, 
and are operated at temperatures between 8 and 10 mK [44,45]. According to the Debye Law, the specific 
heat of Te0 2 crystals is given by C(T) = (3{T/<d D f, where = 1994JK~ 1 mor 1 and 9^ is the Debye 
temperature. In these materials, C{T) is due almost exclusively to lattice degrees of freedom. A special 
measurement determined the value of 9_d, as 232 K [43]. This differs from the previously published value 
of 272 K [46]. The specific heat followed the Debye Law down to 60 mK. The heat capacity of these crystals, 
extrapolated to 10 mK, is 2.3 x 10~ 9 JK _1 . With these values of the parameters, an energy deposition of a 
few keV will result in a measurable temperature increase, AT. In CUORICINO, AT is measured by high- 
resistance germanium thermistors glued to each crystal. More details can be found in reference [44] and in 
earlier publications [47,48]. Accordingly, the temperature increase caused by the deposition of energy equal 
to the total /3/3-decay energy, Q p0 = 2530.3 ± 2.0 keV, would be 1.77 x 10~ 4 K. To obtain usable signals 
for such small temperature changes, very sensitive thermistors are required. 

The thermistors are heavily doped high-resistance germanium semiconductors with an impurity concen- 
tration slightly below the metal-insulator transition. High quality thermistors require a very homogeneous 
doping concentration. CUORICINO uses Neutron Transmutation Doped (NTD) germanium thermistors. 
This is achieved by means of uniform thermal neutron irradiation throughout the entire semiconductor vol- 
ume, in a nuclear reactor. The electrical conductivity of these devices, which is due to variable range hopping 
(VHR) of the electrons, depends very sensitively on the temperature. The resistivity varies with temperature 
according to p = po exp (-^) 7 , where p and T depend on the doping concentration and 7 = 1/2. 

Thermistors can be parameterized by their sensitivity, A(T), defined as follows: A(T) = \d(\nR) /d(lnT)| 
= 7 (Xb/T) 7 , and where the resistance is R(T) = i? exp (T /T) 7 . The parameter R = p (d/a), where d 
and a are the distance between the contacts and the cross section of the thermistor, respectively. The values 
of Rq , To and 7 were experimentally measured for about one third of the thermistors, and the average values 
used for the rest. The measurements were done by coupling the thermistor to a low-temperature heat sink 
with a high-heat-conductivity varnish glue, which can be easily removed with alcohol. The base temperature 
of the heat sink is between 15 and 50 mK [50]. A current flows through the device and an I-V load curve 
is plotted. The curve becomes very non-linear due to the power dissipation, which causes the dynamic 
resistance, the slope of the I-V curve, to invert from positive to negative. The characterization, as discussed 
in Ref. [51] is done on the thermistors directly mounted on a heat sink, while the optimum bias is studied 
for the complete detector, thermistor and crystal, since the noise figure depends on all thermal conductances, 
glue, wires, Teflon, etc. This allows the maximization of the signal to noise ratio. The parameters of each 
thermistor are determined from a combined fit to a set of load curves measured at different base temperatures. 
A detailed description of the characterization process for Si thermistors was described in Ref. [5 1 ] and same 
process was used for the CUORICINO Ge thermistors. 

The thermistors used in the MIBETA and CUORICINO experiments were specially developed and produced 
for this application [52]. It is necessary to optimize the neutron doping of the Ge. This is facilitated by 
foils of metal with long-lived (n, 7) radioactive daughter nuclides, allowing the neutron exposure to be 
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evaluated without having to wait for the intense radiation of the Ge in the Ge sample to decay. Following 
the decay period, the Ge is heat treated to repair the crystal structure and then cut into 3x3x1 mm 
strips. Electrical connections are made with two 50 /im gold wires, ball bonded to metalized surfaces on the 
thermistor. The thermistors are glued to each bolometer by nine spots of epoxy, deposited by an array of 
pins for better control of the thermal conductances and to minimize stresses at the interface between the two 
materials. 



IV THE CUORICINO DETECTOR 



CUORICINO is a pilot experiment for a larger experiment, CUORE (Cryogenic Underground Observatory 
for Rare Events) discussed later. It is a tower of 13 planes [25,26]. As shown in Fig. [T] the CUORICINO 
structure is as follows: each of the upper 10 planes and the lowest one consists of four 5x5x5 cm 3 TeOi 
crystals (of natural isotopic abundance of 130 Te) as shown in the upper right hand figure, while the \\ th and 
12" 1 planes have nine, 3x3x6 cm 3 crystals, as shown in the lower right hand figure. In the 3x3x6 cm 3 
planes the central crystal is fully surrounded by the nearest neighbors for greater veto capability. 







I 


! I.J 


1 



Figure 1: (Color online) The Tower of CUORICINO and individual 4 and 9 detector modules. 

The smaller crystals are of natural isotopic abundance except for four. Two of them are enriched to 82.3% 
in 128 Te and two are enriched to 75% in 130 Te. All crystals were grown with pre-tested low radioactivity 
material by the Shanghai Institute of Ceramics and shipped to Italy by sea to minimize the activation by 
cosmic ray interactions. They were lapped with specially selected low contamination polishing compound. 
All these operations, as well as the mounting of the tower, were carried out in a nitrogen atmosphere glove 
box in a clean room. The mechanical structure is made of oxygen-free high-conductivity copper and Teflon, 
and both were previously tested to be sure that radioactive contaminations were minimal and consistent with 
the required detector sensitivity. 



7 



Thermal pulses are measured with NTD Ge thermistors thermally coupled to each crystal. The thermistors 
are biased through two high-impedance load resistors at room temperature, with resistances typically in ex- 
cess of one hundred times that of the thermistors. The large ratio of the resistances of the load resistors over 
those of the thermistors allows the parallel noise to be kept at an adequate level. Low frequency load-resistor 
noise was minimized by a specially designed circuit [53]. The voltage signals from the thermistors are 
amplified and filtered before being fed to an analog-to-digital converter (ADC). This part of the electronic 
system is DC coupled, and only low-pass anti-aliasing filters are used to reduce the high-frequency noise. 
The typical bandwidth is approximately 10 Hz, with signal rise and decay times of order 30 and 500 ms, re- 
spectively. This entire chain of electronics makes a negligible contribution to the detector energy resolution. 
More details of the design and features of the electronic system are found in [54]. The gain of each bolome- 
ter is stabilized by means of a Si resistor of 50-100 kfi, attached to each bolometer that acts as a heater. 
Heat pulses are periodically supplied by a calibrated ultra-stable pulser [55]. This sends a calibrated voltage 
pulse to the Si resistor. This pulse has a time duration very much shorter than the typical thermal response 
of the detector [44]. The Joule dissipation from the Si resistor produces heat pulses in the crystal almost 
indistinguishable in characteristic shape from those from calibration 7-rays. The heater pulses are produced 
with a frequency of about one in every 300 s in each of the CUORICINO bolometers. Any variation in the 
voltage amplitude recorded from the heater pulses indicates that the gain of that bolometer has changed. The 
heater pulses are used to measure (and later correct offline) for the gain drifts. Two other pulses, one at lower 
and one at higher energies, are sent to the same resistors with much lower frequency. The former is used to 
monitor threshold stability, and the latter to check the effectiveness of the gain stability correction. 

The tower is mechanically decoupled from the cryostat to avoid heating due to vibrations. The tower is 
connected through a 25 mm copper bar to a steel spring fixed to the 50 mK plate of the refrigerator. The 
temperature stabilization of the tower is made by means of a thermistor and a heater glued to it. An electronic 
channel is used for a feed back system [56]. The entire setup is shielded with two layers of lead of 10 cm 
minimum thickness each. The outer layer is made of common low radioactivity lead, while the inner layer 
is made of special lead with a measured activity of 16 ± 4 Bq/kg from 210 Pb. The electrolytic copper 
of the refrigerator thermal shields provides an additional shield with a minimum thickness of 2 cm. An 
external 10 cm layer of borated polyethylene was installed to reduce the background due to environmental 
neutrons. 

The detector is shielded against the intrinsic radioactive contamination of the dilution unit materials by an 
internal layer of 10 cm of Roman lead { 210 Pb activity < 4 mBq/kg [50]), located inside of the cryostat 
immediately above the tower of the array. The background from the activity in the lateral thermal shields of 
the dilution refrigerator is reduced by a lateral internal shield of Roman lead that is 1 .2 cm thick. The refrig- 
erator is surrounded by a Plexiglas anti-radon box flushed with clean iV2 from a liquid nitrogen evaporator 
and is also enclosed in a Faraday cage to eliminate electromagnetic interference. A sketch of the assembly 
is shown in Fig. |2] 

When cooled to 8 mK there is a temperature spread of ~ 1 mK among the different detectors. Routine 
calibrations are performed using two wires of thoriated tungsten inserted inside the external lead shield in 
immediate contact with the outer vacuum chamber (OVC) of the dilution refrigerator. Calibrations normally 
last one to two days, and are performed at the beginning and end of each run, which lasts for approximately 
two-three weeks. 

The CUORICINO array was first cooled down at the beginning of 2003. However, during this operation 
electrical connections were lost to 12 of the 44 detectors of 5 x 5 x 5 cm 3 , and to one of the 3x3x6 cm 3 
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Figure 2: (Color online) A sketch of the CUORICINO assembly showing the tower hanging from the mixing 
chamber and the various heat shields and the external shielding. 

crystals. Thermal stresses broke the electrical connections on their thermalizer stages that allow the transition 
in temperature of the electric signals in several steps from the detectors at ~ 8 mK to room temperature. 
When the cause of the disconnection was found, new thermalizer stages were fabricated and tested at low 
temperature. However, since the performance of the remaining detectors was normal, and their total mass 
was ~ 30 kg, warming of the array and rewiring were postponed for several months while 0^/3/3-decay data 
were collected. At the end of 2003, CUORICINO data acquisition was stopped and the system was warmed 
to room temperature and the broken thermalizer stages were replaced with new ones. During this operation, 
the tower was kept enclosed in its copper box to prevent possible recontamination of the detectors. As a 
consequence, two detectors whose disconnections were inside the box were not recovered. The same was 
true for one of the small central detectors whose Si resistor was electrically disconnected inside the box. 
In the middle of 2004, CUORICINO was cooled down and data collection began again. Typical calibration 
spectra are shown in Fig. [3] 



V DATA ACQUISITION AND ANALYSIS 

The signals coming from each bolometer are amplified and filtered with a six-pole Bessel low-pass filter 
and fed to a 16-bit ADC. The signal is digitized with a sampling time of 8 ms, and a circular buffer is 
filled. With each trigger pulse, a set of 512 samples is recorded to disk; accordingly, the entire pulse shape 
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Figure 3: Typical calibration spectra of the CUORICINO array with a 232 Th source: 5x5x5 cm 3 crystals 
upper frame, 3x3x6 cm 3 crystals lower frame. 



is stored for offline analysis. Each channel (bolometer) has a ompletely independent trigger and trigger 
threshold, optimized according to the bolometers typical noise and pulse shape. Starting with run No. 2, the 
CUORICINO data acquisition (DAQ) now has a software trigger that implements a "debounce" algorithm 
to reduce spurious fast signal triggering. The trigger is ready again within a few tens of ms, a delay due to 
the debounce time. Therefore, most of the pile-up events are re-triggered. The trigger efficiency above 100 
keV was evaluated as 99 ± 1% by checking the fraction of recorded pulser signals. The offline analysis uses 
an Optimal Filter technique [44] to evaluate the pulse amplitudes and to compare pulse-shapes with detector 
response function. Events not caused by interactions in the crystals are recognized and rejected on the basis 
of this comparison. Pile-up pulses are identified and dealt with. This is important for calibration and high 
rate measurements because the pulses have long time durations and pile-up pulses can significantly increase 
the dead time. However, the pile-up fraction during the search for 0^/3/3-decay is negligible given the low 
trigger rate from signals above threshold. The pile-up probability on the rise time is ~ 0.01%, while that on 
the entire sampling window is quite a bit higher, ~ 0.4%. However, these events are easily identified and 
the pile-up pulses are rejected. The total trigger rate, before any pulse-shape rejection, is time and channel 
dependent. On a single channel it ranges from a few mHz to hundreds of mHz, with a mean value of about 
20 mHz. Accepted-pulse amplitudes are then corrected using the variation in the gain measured with the 
heat pulses from the Si resistors. Finally, spectra are produced for each detector. 

Any type of coincidence cut can be applied to the data written to disk, before the creation of the final spectra, 
depending on the specific analysis desired. In the case of /3/3-decay analyses, anticoincidence spectra are 
used. This allows the rejection of background counts from gamma rays that Compton scatter in more than 
one bolometer, for example. The probability of accidental coincidences over the entire detector is negligible 
(< 0.6%). Crosstalk pulses have been observed between a few channels; however, the resulting pulses are 
rejected on the basis of pulse-shape. 
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VI SOURCE CALIBRATION AND DETECTOR 
PERFORMANCE 



The performance of each detector is periodically checked during the routine calibration with the 232 Th 
gamma rays from thoriated calibration wires. The most intense gamma ray peaks visible in the calibration 
spectra are used. They are the: 511, 583, 911, 968, 1588 and 2615 keV 7-rays, and the single escape peak 
of the 2615 keV gamma ray at 2104 keV. The resulting amplitude-energy relationship is obtained from the 
calibration data, and the pulse amplitudes are converted into energies. The dependence of the amplitude on 
energy is fit with a second order log-polynomial for which the parameters were obtained from the calibration 
data. The selection of the functional form was established by means of simulation studies based on a thermal 
model of the detectors. The details of how the thermal model was applied have been published elsewhere 
[44]. These calibration data are also used to determine the energy resolution of each bolometer. Data sets are 
collected for two to six weeks, separated by radioactive-source calibrations. The data collected by a single 
detector in this short time does not have the statistical significance to show the background gamma-ray lines 
because of the very low counting rates. The energy resolution, and the stability of the energy calibration, 
relies on the heater pulses, and on the initial and final source calibration measurements. 

Double-beta decay data collected with each detector during a single data collection period are rejected if any 
of the following criteria are not fulfilled: 

(i) The position of the 2615 keV background 7-ray line from the decay of 208 Tl, in the initial and the 
final source-calibration measurements must be stable to within 1/3 of the measured full width at half 
maximum (FWHM) of the 2615 keV line for that detector. 

(ii) The energy resolution of the 2615 keV 7-ray lines in the initial and final energy calibration measure- 
ments must be stable within 30%. 

(iii) The energy position of the heater pulses during the entire data collection period for that data set must 
be stable to within 1/3 of the characteristic (FWHM) for that detector. 

(iv) The energy resolution measured with the heater pulses for that entire data collection period must be 
stable within 30% over the entire data collection period. 

Whenever one or more of these criteria is not fulfilled, the data from that detector are not included in the 
final data set. Approximately 17% of the data were discarded because they failed one or more of these 
criteria. Frequent causes of failure to satisfy all of the criteria were noise pulses that degrade the energy 
resolution and temperature drifts that change the operating parameters of the bolometers. The particular 
bolometers involved cary; however, some are more sensitive to noise and temperature changes than others. 
The application of coincidence cuts does not change the efficiency; however, the difference in rise time 
between pulses from various bolometers can cause coincidences not to be recognized as such, cut this effect 
is small, in any case, the only result of the failure to recognize coincidences is the loss of background 
reduction, which would tend to make the quoted bound conservative. 

In both runs, the measured detector performances appear to be excellent; the average FWHM resolutions in 
the energy region around 2530 keV during the calibration measurements are 7 and 9 keV, for the 5x5x5 cm 3 
and 3x3x6 cm 3 detectors, respectively. The spread in the FWHM is about 2 keV in both cases. The smaller 
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detectors have somewhat worse resolution on average, while they also exhibit a very important nonlinearity. 
When the calibration spectra from all of the larger and smaller detectors are summed together, the summed 
spectrum resembled that of a single large detector as shown in Figure [3] 




Figure 4: The sum spectrum of the background from the 5x5x5 cm 3 detectors, from both runs, to search 
for 0^/3/3-decay. 



VII DOUBLE-BETA DECAY RESULTS 



Following the shutdown discussed earlier, and restart in May 2004, a second interruption was required 
to remove the malfunctioning helium liquefier used to automatically refill the main bath of the dilution 
refrigerator. There were also short interruptions for routine maintenance of the 17-year old refrigerator. 
Excluding these interruptions, the duty cycle was very satisfactory, ~ 60% , not withstanding the fact that 
15 to 20% of the live time is necessary for calibration. 

The three spectra corresponding to large (5x5x5 cm 3 ) detectors and the smaller natural and enriched 
(3x3x6 cm 3 ) detectors are kept separate because of the different detection efficiencies for /3/5-decay 
events, and also because of their different background counting rates. For similar reasons, the spectra of the 
two runs are treated separately. Because the background rates in the spectra of Runs I and II do not show 
any statistically significant difference, it was concluded that no recontamination of the detector took place 
when the cryostat was opened to air during the interruption between Runs I and II. The full data set used in 
this analysis has a total effective exposure of 11.83 kg • yr of 130 Te for the entire array. 

The full summed spectrum, shown in Fig. |4j clearly exhibits the 7-ray line from the decay of 40 K, and those 
from the 238 U and 232 Th chains. Also visible are the lines of 57 Co, 60 Co, and 54 Mn, due to the cosmogenic 
activation of the tellurium and the copper frame. The correct positions and widths of the peaks in the sum 
spectrum demonstrate the effectiveness of the calibration and linearity of the spectra. The accuracy of 
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calibration in the O^/3/3-decay region was evaluated to be about ±0.4 keV. The details of the gamma-ray 
background resulting from a preliminary analysis of Run 2 are given in Tables ITfl Hill and ITVl There is also 
clear evidence of alpha backgrounds at energies above the 2614.5 keV gamma ray in the decay of 232 Tl. A 
detailed analysis attributes the dominant background in the region of interest to degraded alpha particles on 
the surface of the copper frames. A major effort is underway to reduce this to a minimum. 

Table II: Gamma rays from the decay of 232 Th observed in Run II. 



ijiicrgy ^Kc v ) 




Pr\nntc/1 OHO 
I^OUIILS/ 1WWW 1 


238.6^ 


>m Pb 


6.84 ±0.43 


338." 


228 Ac 


0.89 ±0.40 


463.0* 


228 Ac 


1.33 ±0.25 


510.7^ 


2Q8 Tl 


7.78 ±0.38 


583. ~ 


208 Tl 


3.88 ±0.30 


727.3 


212 Bi 


1.04 ±0.21 


785.4f 


212 Bi 


1.02 ±0.20 


794." 


228 Ac 


0.70 ±0.25 


833.0^ 


228 Ac 


2.85 ±0.25 


911." 


228 Ac 


4.69 ±0.26 


964.8 


228 Ac 


1.37±0.19 


968.9 


228 Ac 


2.79 ±0.21 


1588.1 


228 Ac 


0.65 ±0.12 


1593.0f 


208 Tl 


0.25 ±0.10 


1620.6" 


212 Bi 


0.58 ±0.15 


1631.0 


228 Ac 


0.39 ±0.13 


2614.5 


208 Tl 


6.90 ±0.26 



"Contains a contribution from the U chain. 

''Contains a contribution from 125 Sb. 

'Contains a contribution from annihilation radiation. 

''Contains a contribution from 214 Bi in the U chain. 

'Contains a contribution from 54 Mn. 

-'Contains a contribution from 214 Bi in the U chain. 

The average background counting rates in the region of 0v(3/3 decay are: 0.18 ±0.01, and 0.20 ±0.04 counts 
per keV, per kg, per year (keV^kg^y -1 ) for the 5x5x5 cm 3 and 3x3x6 cm 3 crystals, respectively. The 
sum background spectrum from about 2300 to 2700 keV, of the 5x5x5 cm 3 and 3x3x6 cm 3 crystals, is 
shown in Fig. The shape of the background in the region of interest does not change when anticoincidence 
requirement is applied. An extensive analysis of the background contributions implies that the continuum 
background in the region of interest around 2530 keV breaks down as follows: 10 ± 5% is due to surface 
contamination of the TeOi crystals with 238 U and 232 Th; 50 ± 20% is due to surface contamination of 
the copper surfaces facing the bolometers also with 232 Th and 238 U; and 30 ± 10% is due to the tail of the 
2614.5 keV gamma ray in the decay of 232 Th from the contamination of the cryostat copper shields. Finally, 
there were no observable gamma-ray lines associated with neutron interactions. Monte-Carlo simulations of 
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Table III: Gamma rays from the 238 U chain in data of Run II. Most of the activity is attributed o a radon 
contamination due to temporary leak in the anti-radon box surrounding the refrigerator. 



Energy (keV) 


Isotope 


Rate Cts/1000 h 


Energy (keV) 


Isotope 


Rate Cts/1000 h 


241.91 


2Up b 


6.84 ±0.43 


1401.7 


2U Bi 


1.23 ±0.13 


295." 


214p b 


2.69 ±0.48 


1408.0 


21i Bi 


1.85 ±0.15 


352.0 


214 nl 

Pb 


3.88 ± 0.42 


1509.5 


214 d • 
Bl 


■toe i n n 

1.85 ± 0.13 


609.4 


214 Bi 


13.09 ±0.47 


1583.2 


21i Bi 


0.99 ±0.15 


665.6 


214 Bi 


2.54 ±0.33 


1594.7* 


21i Bi 


0.25 ±0.10 


768.4 


21i Bi 


2.55 ±0.33 


1599." 


21i Bi 


0.43 ±0.90 


786.0^ 


21i Bi 


1.02 ±0.20 


1661.5 


21i Bi 


1.06 ±0.13 


803.0 


2Wp 


1.52 ±0.19 


1729.9 


21i Bi 


2.51 ±0.14 


934.1 


21i Bi 


1.75 ±0.17 


1764.7 


21i Bi 


14.28 ±0.38 


1120.4 


21i Bi 


10.84 ±0.40 


1838.4 


21i Bi 


0.40 ±0.07 


1155.3 


21A Bi 


1.38 ±0.14 


1847.7 


21A Bi 


1.98 ±0.17 


1238.2 


21i Bi 


4.83 ±0.21 


2118.9 


21i Bi 


1.21 ± 0.12 


1281.1 


21i Bi 


1.32 ±0.13 


2204.5 


21i Bi 


4.55 ±0.24 


1377.8 


21i Bi 


3.37±0.17 


2448.0 


21i Bi 


1.51 ±0.14 


1385.3 


2U Bi 


0.88 ±0.88 









"Contains a contribution from Pb in the Th chain. 
^Contains a contribution from 208 Tl in the Th chain. 
' Contains a contributions from 214 Bi in the Th chain. 



Table IV: Background gamma rays from a variety of sources including isotopes produced by cosmogenic 
neutrons: 60 Co, 54 Mn, and fall out isotopes 137 Cs, 207 Bi. 



Energy (keV) Isotope Rate Cts/1000 h Energy (keV) Isotope Rate Cts/1000 h 



T22J ^Co 5.39 ±0.44 661J ^ 1.26 ±0.19 

427.9 125 Sb 1.95 ±0.27 834.8^ 54 Mn 2.86 ±0.25 

463.2* 125 Sb 1.33 ±0.25 1063." 207 Bi 2.36 ±0.29 

511. annihilation 7.78 ± 0.38 1173.2 60 Co 11.6 ±0.33 

569.7 207 B 3.11 ±0.27 1332.5 m Co 11.9 ±0.36 

600.6 125 Sb 1.42 ± 0.20 1 461.0 40 K 31.4 ±0.58 

635.9 125 Sb 0.64 ± 0.18 2 5 05.7 60 Co 0.31 ± 0.05 



"Contains a contribution from 228 Ac in the Th chain. 
^Contains a contribution from 228 Ac in the Th chain. 
c Contains a contribution from 2C18 T£ in the Th chain. 



14 



the neutron shield imply that the background from neutron interactions would be negligible. 

The energy resolution for the complete data set was computed from the FWHM of the 2615 keV background 
7-ray line in the decay of 203 TZ at the end of the thorium chain. The results are 8 keV for the forty operating 
5x5x5 cm 3 crystals, and 12 keV for the eighteen 3x3x6 cm 3 crystals. Clearly visible is the peak at about 
2505 keV due the summing of the 1332.50-1173.24 keV 7-ray cascade in the decay of 60 Co. This is 25.46 
keV, i.e., about 7 sigma of the Gaussian energy resolution peak from the 0i//3(3-decay end-point energy of 
130 Te, and could make a negligible contribution to the region under the expected O^/3/3-decay peak . The 
sum spectrum from 2290 to 2700 keV is shown in Fig. [5] The sum spectrum from 2470 to 2590 keV is 
shown in Fig. [6] 




energy [fcoV] 

Figure 5: The summed background spectrum in the ~ 400 keV region of interest, which includes the 0vf3(3- 
decay energy 2530.3 ± 2.0 keV. 

The details of the operating conditions and parameters of the two CUORICINO data collection periods 
are given in Table IVl The total usable exposure for Run I + Run II is 11.83 kg • yr of 130 Te. The event 
detection efficiencies were computed with Monte-Carlo simulations; they are 0.863 and 0.845 for the large 
and small crystals, respectively. The loss of efficiency of the bolometers is due to beta particles created near 
the surface that escape part of their energy. From the above exposure data we compute: In 2 x Nl x x t — 
2.809 x 10 25 yr, for the large and In 2 x N$ x e$ x t — 4.584 x 10 24 yr for the small crystals. Here, e is 
the detection efficiency, while Nl and N$ are the numbers of 130 Te nuclei in the large and small detectors, 
respectively. 

The /3/3-decay half-life limit was evaluated using a Bayesian approach. The peaks and continuum in the 
region of the spectrum centered on the /3/3-decay energy were fit using a maximum likelihood analysis 
[57,58]. The likelihood functions of six spectra (the sum spectra of the three types of crystals in the two 
runs) were combined allowing for a different background level for each spectrum, and a different intensity 
of the 2505 keV 60 Co sum peak. Other free parameters are the position of the 60 Co peak and the number 
of counts under a peak at the /3/3-decay energy. The same procedure is used to evaluate the 90% CL limit to 
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Table V: Summary of operating parameters for the two CUORICINO data collection periods. From columns 
1 through 8 are listed: the run number, number of large and small detectors, the active mass of 
130 Te, total run time, the calibration time, the time collecting (3(3 - decay data, the total exposure in 
kg-yr, and the usable exposure in kg-yr after rejection of data not fulfilling the quality requirements. 
The total usable exposure is then 11.83 kg • yr. 



Run# 


Detectors 


Active mass 


Run time 


Calibration 


t-00 


Collected 


Used 




large/small 


[kg 130 Te] 


[d] 


[d] 


[d] 


[kg-yr 130 Te] 


[kg-yr 130 Te] 


1 


29/15 


7.95 


240 


24.5 


55.08 


1.2 


1.06 


2 


40/15 


10.37 


983 


108.5 


415.1 


11.79 


10.77 



the number of counts present in the Oz^/3/3-decay peak. 

Assuming Poisson statistics for the binned data, the fit procedure was formulated in terms of the likelihood 
chi-square analysis as described in the following equation: 

6 

Xl = 2 (Vij - n ^ + n ^ ln ( n i,i/yij)) ' 

i=i 

where j indicates the j th spectrum, riij is the number of events in the i th bin of the j th spectrum and yij is 
the number of events predicted by the fit model. 

Fit parameters were estimated minimizing the \\, while limits were obtained, after proper renormalization, 
considering the Xl distribution in the physical region. The response function for each spectrum is assumed 
to be a sum of symmetric gaussian functions, each having the typical energy resolution of one of the detectors 
summed in that spectrum. The experimental uncertainty in the transition energy is considered by means of 
a quadratic (gaussian) term in the above equation. In the region between 2575 and 2665 keV, assuming a 
flat background, the best fit yields a negative number of counts under the peak (—13.9 ± 8.7). However, the 
resulting upper bound on the number of candidate events in the 0z^/3/3-decay peak is n — 10.7 at 90% C.L. 
These values are normalized to a hypothetical sum spectrum of the entire statistical data set in which each 
of the six spectra are weighted according to the corresponding exposure, geometric efficiency, and isotopic 
abundance. The resulting lower limit on the half-life is computed as: 

T\°; 2 ( 130 Te) > ln 2{N L e L + N s e s } t/n (90% CL) 

= (3.27 x 10 25 /10.7) yr = 3.0 x 10 24 yr. 

The dependence of the value of the limit on systematic uncertainties that arise from the method of analyzing 
the data was investigated in detail. These uncertainties reside in the dead time, energy calibration, Q-value, 
and background spectral shape. The main factor influencing the limit is the uncertainty in the background 
spectral shape. 
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For example, changing the degree of the polynomial used to fit the background in the 0^/3/3-decay region 
from to 2 as well as the selection of the energy window used in the analysis can vary the bound from 2.5 
to 3.3 x 10 24 yr. The quoted 90% CL lower bound was computed using the central value, 2530.3 keV of the 
measured double beta decay energy [49]. There is a small dip in the data centered at ~ 2530 keV as shown 
in Figure|5] This has been treated as a statistical fluctuation. 




Figure 6: The total background spectrum from 2470 to 2590 keV. Clearly visible is the sum peak at 2505.74 
keV due to the sum of the 1173.24 and 1332.50 keV 7-ray cascade in the decay of 60 Co. This 
activity is attributed to the 60 Co in the copper frames generated by cosmic ray neutrons while the 
frames were above ground. The solid lines are the best fit to the region using polynomials of the 
order to 2. The three lines in the region of interest are for bounds (68% and 90%) CL on the 
number of candidate /3/3-decay events. 



VIII NUCLEAR STRUCTURE ISSUES 



There is one theoretical viewpoint that holds that the required model space for 130 Te is still very large for 
reliable shell model calculations and must be severely truncated. Accordingly, the quasiparticle random- 
phase approximation (QRPA) is commonly used [59-78]. The results from these calculations, from author 
to author had, until recently, differed significantly for the same nucleus. In Table IVII only the results from 
references [62, 73] differ significantly from the other 13; they correspond to the largest matrix elements. 
In the QRPA approach, the particle-particle interaction is fixed by a parameter, g pp , which is derived in 
various ways by different authors. Two recent papers by Rodin, et ai, give detailed assessments of the 
uncertainties in QRPA calculations of Oz^/3/3-decay matrix elements, and explain many of the reasons for 
the disagreements between the various authors over the years [60, 61]. The numerical values given in these 
articles were corrected in a later erratum [78]. In Table [VTl we list the values of (m„) corresponding to 
T i/2 ( 13 °Te) > 3.0 x 10 24 yr derived using the calculations of various authors. More details are discussed 
later, including the results from recent shell model calculations. 
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Table VI: Various values of {m v ) corresponding to T°J 2 ( 130 Te) = 3.0 x 10 24 yr. 



Authors/Reference Method (m„) (eV) 



[78] Rodin et al, 2007 using 2z//?/3-decay to fix g pp 0.46 

[62] Staudt et al, 1992 pairing (Bohm) 0.19 

[63] Pantis et al, 1996 no p-n pairing 0.52 

[64] Vogel, 1986 0.47 

[65] Civitarese and Suhoen 2006 0.42 

[66] Tomoda, 1991 0.42 

[67] Barbero, et al, 1999 0.33 

[68] Simkovic, 1999 pn - RQRPA 0.68 

[69] Suhoen et al, 1992 0.64 

[67]Mutoef al, 1989 0.39 

[71] Stoicaef al, 2001 0.60 

[72] Faessler et al, 1998 0.55 

[73] Engel et al, 1989 seniority 0.29 

[74] Aunola et al, 1998 0.41 

[79] Caurier et al, 2008 Nuclear Shell Model 0.58 



Extracting the effective Majorana mass of the electron neutrino from the half-life requires the calculation of 

the nuclear structure factor, fjy = G 0l/ \ Mp' — {$aI 9F*f A/gyj, in Eq. (0. This is not straightforward 

for the nuclei that are the best candidates for 0z//?/3-decay experiments, e.g., 130 Te, because they have many 
valence nucleons. To create a tractable shell-model calculation for these heavy nuclei, it is necessary to 
truncate the model space to the point that could affect the reliability of the results. Accordingly, schematic 
models are employed. As stated above, QRPA has become the standard approach for both 2vf3(3 and Qv(3j3 
decay. The results calculated with QRPA, however, depend on the selection of a number of parameters, and 
the fact that different authors select the parameters in various ways has resulted in large differences in the 
resulting matrix elements as discussed in Ref. [61]. 

In Table|Vl] we list 14 different values of (m v ) derived with QRPA and with renormalized QRPA, (RQRPA), 
corresponding to T®" 2 ( 130 Te) = 3.0 x 10 24 yr, and also the recent shell-model calculations of Caurier et 
al [79]. From the table it is clear that the different ways of applying the same basic model has lead to 
a spread in the resulting matrix elements, and hence in the corresponding value of (rn^), of a factor of 
three [61-74]. This corresponds to differences of a factor of nine in the predicted half-life for a given value 
of (m v ), if all calculations are given the same weight. This assumption, however, cannot be justified. It 
should be recognized that calculation techniques, as well as computational power have made significant 
progress over the years, improving the reliability of both QRPA and shell-model calculations. 

In their recent article, Rodin, Simkovic, Faessler, and Vogel (Tubingen) [61], give detailed discussions of 
how the choices of various parameters in similar models can lead to such discrepancies. These are the gap 
of the pairing interactions, the use of (renormalized) RQRPA that partially accounts for the violation of the 
Pauli principle in the evaluation of the two-fermion commutators, the nucleon-nucleon interaction potential, 
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the strength of the particle-hole interactions of the core polarization, the size of the model space, and the 
strength of the particle-particle interaction, parameterized by the quantity g pp . The matrix elements of the 
virtual transitions through states with J" = 1+ in the intermediate nucleus are extremely sensitive to the 
value of g pp , which makes 2^/3/3-decay matrix elements also very sensitive to it because this decay mode 
only proceeds through 1 + intermediate states. On the other hand, O^/3/3-decay also proceeds via higher 
multipoles through states of higher spin. These transitions are found to be far less sensitive to the value 
of gpp. For this reason, Rodin et al. select the value of g pp that makes the calculation of the 2^/3/3-decay 
half-life agree with the experimental value. In addition, some calculations are greatly simplified by using 
an average energy in the denominator of the second-order matrix-element expression, and the sum over the 
intermediate states is done by closure. When the value, g^ = 1.245, of the axial-vector coupling constant 
obtained from muon decay is used, it commonly lead to a value of the Gamow-Teller strength typically larger 
than the measured value. To ameliorate this situation, a quenched value gA = 1-00 is used. In calculated 
rates of 2^/3/3-decay, which proceed only through J 71 ' = 1 + states, this results in a factor of 2.44 reduction 
in the rate. Using the technique of Rodin et al. [61], the choice of g a A = 1.00 reduces the rate by between 
10 to 30%, depending on the particular nucleus. 

Another serious difference between some of the Of /3/3-decay calculations is due to the treatment of the 
short-range correlations in the nucleon-nucleon interactions. It was also pointed out by Simkovic et al. [68], 
that including the momentum dependent higher order terms of the nucleon current typically result in a 
reduction in the calculated value of the Ov /3/3-decay matrix element by about 30%. These were included in 
the calculations ofRefs. [60,61]. 

In recent paper by Alvarez et al. [75], a QRPA formalism for 2v /3/3-decay in deformed nuclei was presented. 
A considerable reduction in the matrix elements was observed in cases in which there was a significant 
difference in the deformations of the parent and daughter nuclides. Exactly how this would affect 0v(3f3- 
decay calculations is not yet clear. It must be understood that this uncertainty, when resolved could result 
in a further reduction in neutrinoless double-beta decay matrix elements calculated within the framework of 
QRPA and RQRPA. 

In general, however, the paper by Rodin et al. [61], represents a detailed study of the various factors that 
cause the large variations in the nuclear matrix elements of Ov /3/3-decay calculated by different authors over 
the years, and must be taken seriously. The procedure of Rodin et al. [59-61] has the attractive feature 
that it gives a straightforward prescription for selecting the very important particle -particle parameter, g pp . 
However, Civitarese and Suhonen (referred to as the Jyvaskyla group) have given strong arguments in favor 
of using single /3 ± -decay and electron capture data for this purpose, while giving arguments against using 
experimental 2v /3/3-decay half lives [65]. They argue that only states with spin and parity 1+ can be the 
intermediate states involved in 2v /3/3-decay, and that in the neutrinoless process these states play a minor 
role, and that the higher spin states play a dominant role. The Jyvaskyla group recently presented a preprint 
in which they show that the effects of short-range correlations have been significantly overestimated in the 
past [76,77]. Accordingly, their matrix elements originally gave a very different picture of the of the physics 
impact of the CUORICINO data presented in this paper. However, recently there have been some very 
important developments discussed below. 
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IX RECENT DEVELOPMENTS IN 
QRPA CALCULATIONS 



We adopt the position that the large dispersion in values in the nuclear matrix elements implied by the 
values in Table |VT] does not reflect the true state of the art. Instead, we assume that there has been sig- 
nificant progress in understanding the key theoretical issues, as well as large increases in available com- 
putational power. Until very recently, however, two of the recent extensive theoretical treatments of the 
Oz//3/3-decay matrix elements disagreed significantly, and in particular in the case of 130 Te. The relevant 
nuclear structure factors, Fjy, from the Jyvaskyla and Tubingen groups for qa = 1-25 were Fn ( 130 Te) = 
1.20 ± 0.27 x 10~ 13 yr" 1 of Rodin et al. [61], and F N ( 130 Te) = 5.13 x 10~ 13 yr" 1 of Civitarese and 
Suhonen [65]. 

Recently an erratum was submitted by Rodin et al. [78] with major corrections to Table U of Ref. [61]. A 
coding error was discovered in the computation of the short-range correlations that, for example, increased 
the predicted Of /?/3-decay rate of 130 Te by a factor of 4.03. Their corrected value of the nuclear structure 
factor of 130 Te, is now Fn ( 130 Te) = 4.84^^4 x 10~ 13 yr -1 , in good agreement with the above value 
given by Civitarese and Suhonen. However, there is still a small disagreement between these two groups 
concerning the technique for calculating short-range correlations. Rodin et al, used a Jastrow-correlation 
function, which has subsequently been shown by Kortelainen et al. [76] to overestimate the effects of short- 
range correlations, and hence to result in an excessive reduction in the nuclear matrix elements. 

Kortelainen et al. [77] have also updated the calculations of Civitarese and Suhonen. They extended their 
model space, for the cases of 116 Cd, 12 M30 j> e an( j ^Xe, to include the lp-0f-2s-ld-0g-2p-lf-0h single 
particle orbitals, calculated with a spherical Coulomb-corrected Woods-Saxon potential. In Ref. [77], a 
complete discussion is given of their method of fixing the parameters of the Hamiltonian. In this treatment 
they fix particle-particle parameter g pp of the pnQRPA using the method of Rodin et al. [59-61], namely 
with the experimentally measured 2z//3/3-decay half-lives. They did not use the Jastrow-correlation function 
to correct for the short-range correlations, but rather they employ a "unitary correlation operator method" 
(UCOM), which in the case of 130 Te increases the matrix element by a factor of 1.38 over that calculated 
with the Jastrow correlation function. Their new values for the nuclear structure factors are: 



F N ( 130 Te) , or = 7.47 x 10~ 13 yr" 1 , 
^( 130 T e ) M=LOO =4.93xlO- 13 yr- 1 . 

This is to be compared to the results of the earlier work of Civitarese and Suhonen [65]. 

In any case, the major disagreements between the Jyvaskyla and Tubingen groups have finally been under- 
stood, and the present difference in the predicted 0^/3,3-decay rates of 130 Te now differ by a factor of 1.06, 
whereas the earlier disagreement was by a factor of 4.28. Some remaining differences might well lie in the 
differing methods of applying the short-range correlations (see also the discussion in Ref. [80]). In any case 
these recent developments have had a major impact on the interpretation of the CUORICINO data. 

Furthermore, the group of Caurier et al. [79], have recently given new values for these matrix elements from 
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improved nuclear shell model calculations. The shell-model matrix elements are somewhat smaller than 
those of the recent Jyvaskyla and corrected Tubingen results, and according to their matrix elements, the 
CUORICINO data imply: (m„) < 0.58 eV. 



The CUORICINO array is the only operating 0;//5/3-decay experiment, with energy resolution adequate 
to potentially probe the range of effective Majorana mass, (m u ), implied by the observation of Qvj3j3- 
decay claimed by Klapdor-Kleingrothaus et al. [19,20]. In the 2006 article by Klapdor-Kleingrothaus and 
Rrivosheina (KK& K) [20], the peak in the spectrum centered at Qpp = 2039 keV is interpreted as the 
0i/^/3-decay of 76 Ge, consistent with the range: T°J 2 ( 76 Ge) = {1.30 - 3.55} x 10 25 yr (3a). The best-fit 

value is (2.231q' 31 ) x 10 25 yr. In this discussion we offer no critique of the claim, however, since this claim 
has been criticized from several points of view [21-23], it is interesting to ask if it is feasible to observe 
a 0i//3/3-decay with this half-life with a significant confidence level with the published parameters of the 
experiment. Below, we show that the answer is "yes", the experiment could have made the observation in 
the range of half-lives quoted [20]. 

It is straightforward to derive an approximate analytical expression for the half-life sensitivity for discovery 
at a given confidence level that an experiment can achieve (see Appendix). The achievable discovery half- 
life, when the background rate is nonzero, is expressed as: 



It is more conventional to simply have b 6(E) in the denominator of the root of Eq. (O as prescribed by the 
Particle Data Book [81]. However, when the background continuum is obtained by a best fit to all peaks 
and continuum in the region, we choose this alternative approach. In Eq. (O, n a is the desired number of 
standard deviations of the (CL) (3 for CL = 99.73%, for example), e is the event detection and identification 
efficiency, a is the isotopic abundance, W is the molecular weight of the source material, M is the total mass 
of the source, ( is the signal-to-background ratio, b, is the specific background rate in counts/ke V/kg/yr, and 
6(E) is the instrumental width of the region of interest related to the energy resolution at the energy of the 
expected Oz^/3/3-decay peak. 

The values for these parameters for the Heidelberg-Moscow experiment [17, 19, 20] are: Mt = 71.7 kg • yr, 
b = O.llkg^keV^yr- 1 , e = 0.95, a = 0.86, W = 76, and 6(E) = 3.27 keV. The number of counts 
under the identified peak at 2039 keV is 28.75 ± 6.86. The average value of the background near the region 
of interest was 11.6 counts, therefore £ = 2. Direct substitution into Eq. (0 yields: 



X CUORICINO AS A TEST OF THE CLAIM 
OF DISCOVERY 




(9) 



T i/2 76 Ge) = 0.9 x 10 25 yr; Ty 2 (3a) = 1.2 x 10 25 yr. 



(10) 
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Using the less conservative approach with bS(E) in the denominator, the predicted half-life sensitivity for a 
discovery is 



(4a, 



76 



Ge) = 1.6 x 10 25 yr; 



(3<r) 



2.13 x 10 25 yr. 



(11) 



These are close to the claimed most probable value given in Ref. [20]. This simple analysis is independent 
of the claimed result, with the exception of the determination of the signal to background ratio, £. The con- 
clusion is that with the given experimental parameters, this experiment could have had a discovery potential. 
Since this analysis does not account for statistical fluctuations, the discovery confidence level could possibly 
fall between 3cr and 5cr. Any criticism of the claim would involve a reanalysis of the data, and the interpre- 
tation of the background peaks in the region. This falls outside of the scope of this discussion. Accordingly, 
we do not question the claim, but rather ask how well the present CUORICINO data confront it, now and in 
the future after five years of running. 

While the many theoretical calculations of the nuclear matrix elements over the years have differed signifi- 
cantly, the recently corrected-QRPA calculations of Tubingen [78], those of Jyvdskyla [65], and shell model 
calculations of Caurier et al. [79], differ by less than about 30%. We have chosen to use for further analysis 
of the physics impact of the present CUORICINO data. 

Equation ([H) can be inverted to obtain the values of the nuclear structure factor, Fn, using the calculated half- 
lives for Oz^/3/3-decay calculated with a given (m u ) by the authors of the theoretical papers. The resulting 
values are as follows: 



76 Ge 



fl A =1.245 • 



Rodin, et al. : 
Caurier, et al.: 
Civitarese and Suhonen: 



F N = 1.22 



+ 0.10 



x 10 



-13 , 



F N = 4.29 x 10 
F N = 7.01 x 10~ 14 



(12) 



yr 



S 3 A = 1.245- 

Rodin, et al. : 
Caurier, et al: 
Civitarese and Suhonen: 



4.84±i-j° x 10" 



Fn 

F N = 2.57 x 10" 13 
F N = 5.13 x 10 



13 



yr 1 (corrected value), 



yr 



13 yr" 1 . 



(13) 



The resulting values and ranges of values of (m„) implied by the KK&K data, and by the CUORICINO data 
are as follows: 
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<^>kk & k = {0-23- 0.43} eV, 
Mcuor < {0.38 - 0.46} eV, 

(m,)^ &K = {0.32 - 0.54} eV, (14) 

(«V)cuOR ^ °' 41 eV ' 

<^>kk & k = {0-41 - 0.68} eV, 

("V>CUOR < °- 58 eV - 

The results of the analyses with the new corrected matrix elements of Ref. [78] imply that the CUORICINO 
sensitivity has entered well into the range of values of (m„) implied by the claim of KK&K. In the other 
two analyses, the CUORICINO data also constrain part of the range of values of (m v ) implied by KK&K. 



It is also interesting to try to predict the sensitivity of CUORICINO if it were to continue to operate for a 
total of 5 years. The three recent calculations of the nuclear matrix elements result in the following predicted 
decay rates if the Heidelberg claim is correct. In this case, the decay rates would be: 



^kk&k (j Ge) = {1.95 - 5.32} x lO" 26 yr" 1 , 
T Roa ( 13 °Te) = {0.62 - 2.94} x 10~ 25 yr" 1 , (15) 
T civ ( 13 °Te) = {1.43 - 3.89} x 10~ 25 yr" 1 , 
T cl ( 13 °Te) = {1.17 - 3.19} x 10~ 25 yr" 1 . 

Accordingly, we can calculate the number of O^/3/3-decay counts with 5 years of live-time operation expected 
in the CUORICINO data consistent with the claim of KK&K. The exposure would be: Nte = 2.85 x 10 26 
y, resulting in the following predicted number of real O^/3/3-decay events: 



{18-84V/3/3, 

{41 - 110} , w , (16) 
{33 - 91V/3/3. 

These counts would be superimposed on an expected background of 35 to 39 counts per keV in the 8 keV 
region of interest centered at 2530 keV. 

The constraints placed by the current CUORICINO data might favor the lower numbers in the ranges above. 
This would make it more challenging for CUORICINO to confirm the discovery claim of KK&K, and 
renders it almost impossible to rule out the KK&K claim with a significant level of confidence. The solution 
to this problem is the construction and operation of the proposed first tower of CUORE, called CUORE-0, 
combine its data with that of CUORICINO, and later the complete CUORE Experiment. 



r^Nte 



'Cau 



Nte 
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XI THE PROPOSED CUORE EXPERIMENT 



The proposed CUORE detector will be made of 19 towers of TeC?2 bolometers, very similar to the CUORI- 
CINO tower [28]. Each will house 13 modules of four 5x5x5 cm 3 crystals with masses of ~ 750 g. 
CUORE will contain — 200 kg of 130 Te. The 988 bolometers will have a total detector mass of ~ 750 
kg and will operate at 8-10 mK. An intense research and development program is underway to reduce the 
background to 0.01 counts/keV/kg/yr. Thus far a reduction has been achieved that has reached within a 
factor of 2.4 of this goal in the region of 2530 keV, the Q-value for the 0^/3/?-decay of 130 Te. With this 
background, CUORE would reach a sensitivity of ~ T-P^ ( 130 Te) w 2.1 x 10 26 yrin 5 years. The secondary 
goal is to achieve a background level of 0.001 counts/keV/kg/yr. This would allow a half-life sensitivity of 
T°Y 2 w 6.5 x 10 26 yr. 

In case that the background would be reduced to 0.001 counts/keV/kg/yr, the associated sensitivities in the 
effective Majorana mass of the electron neutrino, (m v ), would be 



(m,) Rod = {0.026- 0.031} eV, 

K) a = 0.028 eV, (17) 

(«^>Cau = °- 040 



The half-life sensitivity is directly proportional to the abundance, a, of the parent /3/3-decay isotope [see 
equation ©]. Accordingly, enriching the detectors of CUORE from 33.8% in 130 Te to 90%, CUORE 
would achieve the same sensitivity with a background of 0.01 counts/keV/kg/yr as it would with natural Te 
and a background of 0.0014 counts/keV/kg/yr. An R&D program, to determine the feasibility and cost of 
isotopically enriching CUORE is underway . In addition, the CUORE collaboration has a rigorous R&D 
program to improve the energy resolution from an average of 8 keV, as it is in CUORICINO, to 5 keV. This 
resolution should be achievable because some of the CUORICINO bolometers have already achieved 5 keV. 
An intense program is underway to determine the cause of the spread in energy resolution. If in the end, 
CUORE does achieve the background of 0.001 counts/keV/kg/yr, in addition is enriched, and has an average 
energy resolution of 5 keV, it could reach a half life sensitivity of 2.5 x 10 27 yr in 10 years. 

In this case the sensitivities become: 



( m ^>Rod = {13- 16}meV, 

(m„) Civ = 14meV, (18) 

( TO ^)cau = 20 meV - 



This brings the sensitivity into the normal hierarchy region, which exceeds the goals of some of the other next 
generation experiments. It is possible to proceed as planned with a natural abundance version of CUORE, 
and then the bolometers could be replaced with those isotopically enriched in 130 Te. This would increase 
the half-life reach by a factor of 2.5 for an enrichment of 85%. 
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XII SUMMARY AND CONCLUSIONS 



The CUORICINO detector is an array of 62 TeC>2 bolometers operating at a temperature of about 8 mK. 
It has a total mass of 40.7 kg of TeO^, containing 11 kg of 130 Te. It has operated for a total exposure 
of N ( 130 Te) te = 5.47 x 10 25 yr, with no observation of 0^/3/3-decay events, results in a lower bound, 
T®f 2 ( 130 Te) > 3.0 x 10 24 yr. The corresponding upper bound on the effective Majorana mass of the 
electron neutrino, (m„), using the corrected nuclear structure calculations of Rodin etal, is (m u ) < (0.38 — 
0.46) eV, while using those of Civitarese and Suhonen yields (m„) < 0.47 eV. With the recent shell model 
calculations the CUORICINO data imply (m„) < 0.58 eV. In all cases, the present CUORICINO data probe 
a significant portion of the range of the half life measured by KK&K. If the Heidelberg claim is correct, the 
nuclear structure calculations of Ref. [78] imply that after 5 years of live time CUORICINO would detect 
{18 — 84}, 0^/3/?-decay events, while those of Ref. [65] imply it would detect {41 — 110} events, and those of 
Ref. [79] imply it would detect {33 — 91} O^/3/3-events. In all cases, these counts would appear in Gaussian 
peaks with FWHM = 8k keV, superimposed on an average background of 35 — 39 counts keV -1 . 

In any case, the current results imply that the continued operation of CUORICINO is very important since 
it represents the only possibility of testing the claim of evidence of O^/3/3-decay for the next 5 years or 
more. 
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A APPENDIX 

An approximate expression for estimating the 0^/3/3-decay half-life at which a given experiment can achieve 
discovery at the confidence level corresponding to n a u, can be derived by reference to Figure [7] Let "C" 
be the total number of counts found in the region of the expected 0^/3/3-decay peak; let "B" be the total 
number of background counts in the same energy interval, 8(E). For the number of real 0^/3/3-decay events 
to have a statistical significance of n a , the following must be true: C — B = n a \[C. In the usual case where 
B^0,a desired signal to background ratio, £ = (C — B)/B, can be chosen; hence C = (1 + £) B. The 
usual expression for the corresponding half-life can be written in terms of these parameters as: 
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Ty 2 (n 



(In 2) ATte 



i a y/(l + B 



(A.l) 



where A~ is the total number of parent nuclei, e is the total detection efficiency, and t is the live time of the 
data collection. The number of parent nuclei can be written in terms of, M, the total mass of the source 
(in an oxide for example), as follows: A" = (lO 3 g/kg/Wq/mole) • (Aq at/mole) ■ a (abundance) • Afkg. 
Substituting these values, and expressing the background in terms of the background rate, B = bM5(E)t, 
where b = (counts/keV/kg/yr), the expression is written: 



t\ ; 2 K) 



4.17 x 10 




bS(E) 



(A.2) 



8(E) 



C'Ti 



B-ii 



Figure 7: Diagram showing the scheme on which Eq. ( IA.21 > is derived 



Of course in the case of zero background, Eq. ( IA.U is used, and the quantity, (1 + £) B, is replaced the 
number of real events in the peak. In case there are no real or background events, i.e., C = B = 0, the 
denominator of Eq. ( IA.U is replaced by the usual quantity, In {1/(1 — CL)}, which is 2.3, (90% C.L.) for 
example, and T®J 2 becomes an experimental lower limit. In Eq. ( IA.2I ). we use the fluctuation in the real 
events instead of that of the background because in these experiments the background level used is that of a 
best fit curve to the background in the region, and the fluctuation is a fitting error and is much smaller than 
the statistical fluctuations in the region of interest. 
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